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The phenomenon of real-space condensation is encountered in a variety of situations such as 
aggregation and fragmentation processes, granular clustering, phase separation, traffic and networks. 
Unlike traditional Bose- Einstein condensation in the momentum space, a condensate in these systems 
forms in real space, e.g., upon increasing the density beyond a critical value a macroscopically large 
mass/cluster may form at a single site on a lattice. In this brief review, I discuss some recent 
, developments in understanding the physical and mathematical mechanism behind this real-space 

' condensation in a class of simple stochastic mass transport models. 
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I. INTRODUCTION 



The phenomenon of Bose-Einstein condensation (BEC) in an ideal Bose gas is by now a textbook material and has 
,_i recently seen a huge revival of interest driven mostly by new experiments. Consider an ideal gas of N bosons in a 
d-dimensional hybercubic box of volume V = L"^. In the thermodynamic limit N —> oo, — > oo but with the density 
O N/V fixed, as one reduces the temperature below a certain critical value Tc{p) in d > 2, macroscopically large 

^ number of particles (oc V) condense on to the ground state, i.e., in the zero momentum quantum state. Alternately, 
G , one encounters the same condensation transition upon fixing the temperature but increasing the density p beyond a 
4!) i critical value Pc{T). 

^ I The traditional BEC happens in the momentum (or equivalently energy) space. In contrast, over the last two 
c/2 . decades it has been realized that a 'similar' Bose-Einstein type condensation can also occur in real space in the steady 
state of a variety of nonequilibrium systems such as in cluster aggregation and fragmentation jamming in traffic 
and granular flow 0, [3| and granular clustering [![ . The common characteristic feature that these systems share is 
the stochastic transport of some conserved scalar quantity which can simply be called mass. Condensation transition 
in these systems occurs when above some critical mass density a single 'condensate' captures a finite fraction of the 
^ total mass of the system. 'Condensate' corresponds to a dominant cluster in the context of granular clustering, or a 
Q ' single large jam in the context of traffic models. Another example of condensation is found in the phase separation 
O . dynamics in one dimensional driven systems where the condensation manifests itself in the emergence of a macroscopic 
domain of one phase [5]. Other examples of such real-space condensation can be found in the socioeconomic contexts: 
I , for example, in wealth condensation in macroeconomics where a single individual or a company (condensate) owns a 
^ ■ finite fraction of the total wealth or in growing networks where a single node or hub (such as 'google') may capture 

' a finite fraction of links in the network 0] ■ 
G^ This real -space condensation mentioned above has been studied theoretically in very simple stochastic mass trans- 

port models defined on lattices. These are typically nonequilibrium models without any Hamiltonian and are defined 
by their microscopic dynamical rules that specify how some scalar quantities such as masses or a certain number of 

■ particles get transported from site to site of the lattice. These rules typically violate detailed balance Under these 
, rules, the system evolves into a stationary or steady state which are typically not Gibbs-Boltzmann state as the system 

■ lacks a Hamiltonian Q. For a class of transport rules, the system can reach a steady state where upon increasing the 
density of mass or particles beyond a critical value, a macroscopically large mass (or number of particles) condenses 
onto a single site of the lattice, signalling the onset of 'real-space' condensation. 

In this article we will mostly focus on homogeneous systems where the transport rules are independent of sites, 
i.e., the system is translationally invariant. In the condensed phase and in an infinite system, the condensate forms 
d ' at a single site which thus breaks the translational invariance spontaneously. In a finite system, the condensate at 
a given site has a finite lifetime beyond which it dissolves and then gets relocated at a different site and the various 
time scales associated with the formation/relocation of the condensate diverge with the increasing system size (see 
later). In heterogeneous systems where the transport rules may differ from site to site, the condensate may form at 
a site with the lowest outgoing mass transport rate [13, El • The mechanism of the condensation transition in 
such heterogeneous systems is exactly analogus to the traditional BEC in momentum space and the site with the 
lowest outgoing mass trasfer rate plays the role of the ground state in the quantum system of ideal Bose gas. In 
contrast, the mechanism of condensation in homogeneous systems, the subject of focus here, is rather different: the 
onset and formation of a condensate in an infinite system is associated with the spontaneous breaking of translational 
invariance. Also, unlike the traditional equilibrium Bose gas in a box, this real-space condensation in nonequilibrium 
mass transport models can occur even in one dimension. 

The purpose of these lectures is to understand the phenomenon of real-space condensation in homogeneous systems 
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within the context of simple one dimensional mass transport models. The main questions we will be addressing are 
threefold: (i) When does the condensation happen- i.e. to find the criterion for condensation (ii) How does the 
condensation happen-i.e., to unfold the mathematical mechanism behind such a transition if it happens and (iii) 
What is the nature of the condensate- e.g., to compute the distribution of mass or the number of particles in the 
condensate. 

The article is organized as follows. In Section II we will discuss three simple and well studied lattice models of 
stochastic mass transport. In Section III, we will consider a generalized mass transport model that includes the 
previous three models as special cases and investigate its steady state. In particular, we will study in detail steady 
states that are factorisable. The necessary and sufficient conditions for such factorisable property will be discussed. 
Thanks to this factorisable property, a detailed analytical study of the condensation is possible for such steady 
states which will be illustrated in Section IV. In Section V we will illustrate how various results associated with 
the condensation transition in factorisable steady states can be simply understood in terms of sums and extremes 
of independent and identically distributed (i.i.d) random variables. Finally we will conclude in Section VI with a 
summary and other possible generalizations/issues associated with the real-space condensation. 



II. THREE SIMPLE MASS TRANSPORT MODELS 



A. Zero Range Process 

The Zero Range Process (ZRP), introduced by Spitzer [T^, is perhaps one of the simplest analytically solvable 
model of mass / particle transport that exhibits a real-space condensation in certain range of its parameters-for a 
review see hj. ZRP is defined on a lattice with periodic boundary conditions. For simplicity, we will consider a 
1-d lattice with L sites, the generalization to higher dimensions is straightforward. On each site of the lattice at any 
instant of time rests a number of particles, say rrii at site i where > is a nonncgative integer. We can also think 
of each particle carrying a unit mass, so that represents the total mass at site i. A configuration of the system 
at any given instant is specified by the masses at all sites {mi, m2, . . . , Wi}. One starts from an arbitrary initial 
condition with a total mass M = m^. The subsequent dynamics conserves this total mass or total particle number, 
or equivalently the density p = M/ L. 

The system evolves via continuous-time stochastic dynamics specified by the following rules: 

• In a small time interval dt, a single particle from site i with number of particles is transported to its right 
neighbour i + 1 with probability U{mi)dt provided rrii > 1. In terms of mass, this means a single unit of mass is 
transferred from site i to site i + 1 with rate U{mi). 

• Nothing happens with probability 1 — U{mi)dt. 

Here U{m) is an arbitrary positive function, with the constraint that U{0) — 0, since there can not be any transfer of 
unit mass if the site has no mass at all. Thus in ZRP, the particle or mass transfer rate t/(m) depends only on the 
number of particles/mass m at the departure site prior to the transfer. One can of course generalize easily the ZRP 
to discrete-time dynamics, with symmetric transfer of particles to both neighbours etc [l^. But here we stick to the 
asymmetric continuous-time model for simplicity. 

As the system evolves under this dynamics, the probability of a configuration P(toi, TO2, ■ • ■ , evolves in 
time and in the long time limit, t ~^ oo, it approaches a time-independent stationary joint distribution of masses 
P(mi, 7712, . . . , mi,). This is the basic quantity of interest, since the statistics of all other physical observables in 
the steady state can, in principle, be computed from this joint distribution. In many such nonequilibrium systems, 
computing the steady state P{mi^m2, ■ ■ ■ ,777/,) is, indeed, the first big hurdle Fortunately, in ZRP, this can be 
computed explicitly and has a rather simple factorised form [T^ . [T3 | 



/(777l)/(7772) . . . /(?77l) S 



P(mi,m2,...,7r7L) = ^^^n,r\ f{'nii)f{m2)...f{mL)5 [} ^rm-M] (1) 
where the weight function /(777) is related to the transfer rate U{m) via 

rn ^ 

f{m) ^ Y[ JJ(j^ foi" m>l 

= 1 for 777 = (2) 
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The delta function in Eq. ([T]) specifies the conserved total mass M and Zl{M) is just a normalization factor that 
ensures that the total probability is unity and satisfies a simple recursion relation 
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Zl{M) = ^ n /(™») M E ™» " = E fi-m^ZL-iiM - m). (3) 

rui i=l \ i / m=0 

To prove the result in Eq. ([T]) one sirnply writes down the Master equation for the evolution of the probability in the 
configuration space and then verifies |1J| that the stationary solution of this Master equation is indeed given by ([T]). 

Finally, the single site mass distribution p{m), defined as the probability that any site has mass m in the steady 
state, is just the marginal obtained from the joint distribution 

t \ s:^ -nl \ J-/ x-^L-l(M-m) 

p(m) ^ > P{m,m2,m3,...,mL) = f{m) „ , . . (4) 

^ ZL[m) 

m2,m3,...,mL 

Note that p{m) implicitly depends on L but this L dependence has been suppressed for notational simplicity. This 
single site mass distribution is important as any signature of the existence of a condensate will definitely show up in 
the explicit form of p{m) . 

Evidently the steady state p{m) depends on the transfer rate U{m) through the weight function /(m) in Eq. ([2]). 
Not all choices of U{m) lead to a steady state with a condensation transition. Indeed, one may ask what choices of 
U{m) may lead to a condensation transition. An example of such a choice is given by U (m) cx (1 + 7/771) for large m, 
which leads to, using Eq. a power law weight function, f{m) ~ m^'' for large m. In this case, it was shown |3, lisj 
that for 7 > 2, the system undergoes a condensation transition as one increases the density p through a critical value 
Pc 1/(7— 2). The condensation transition shows up in p(m) in the thermodynamic limit, which has different forms 
for p < pc, p = Pc and p > Pc [13] 

p{m) ^ exp [—m/m*] for p < pc 

~ for p = Pc 

~ — — [-"condensate" for p > Pc (5) 

Thus for p < Pc, the single site mass distribution decays exponentially with a characteristic mass m* that diverges 
as p Pc from below, has a power law form exactly at p = pc and for p > pc, while the power-law form remains 
unchanged, all the additional mass (p — Pc)L condenses onto a single site which shows up as a bump in p{m) at the 
tail of the power law form (see Fig. ([3])). The term "condensate" in Eq. ^ refers to this additional bump. Physically 
this means that a single condensate coexists with a background critical fluid for p > pc- This change of behavior of 
p(m) as one increases p through pc is a prototype signature of the real-space condensation and one finds this behavior 
in various other stochastic mass transport models that will be discussed below. In addition, many details of the 
condensation phenomena in ZRP also follow as special cases of the more general mass transport model defined in 
Section III. 

Before ending this subsection, it is useful to point out that there have been several other issues and studies on ZRP 
and related models that arc not covered here. The interested readers may consult the reviews [isl [l^. 



B. Symmetric Chipping Model 

Here we discuss another simple one dimensional mass transport model that also exhibits a condensation phase 
transition in its steady state. As in ZRP, this model is also defined on a lattice with periodic boundary conditions 
where each site i carries an integer mass > A non-lattice mean-field version of the model was studied in 
Ref. [1^ . The system evolves via the continuous- time dynamics defined by the following rules [l], [l^ : 

• diffusion and aggregation: in a small time interval dt, the entire mass rrii from site i moves either to its right 
neighbour {i + 1) with probability dt/2 or to its left neighbour {i — 1) with probability dt/2. 

• chipping: in the same interval dt, only one unit of mass chips off site i with mass rrii (provided rrii > 1) to either 
its right neighbour with probability wdt/2 or to its left neighbour with probability wdt/2. 

• With probability 1 — {1 + w) dt nothing happens. 
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Once again the total mass M = pi is conserved by the dynamics. The model thus has two parameters p (density)and 
w (the ratio of the chipping to the diffusion rate). At long times, the system evolves into a steady state where 
the single site mass distribution p{m), for large L, exhibits a condensation phase transition at a critical density \M , 
Pc{w) = \/w + 1 — 1. Remarkably, this equation of state Pc{w) = ^/w + 1 — 1 turns out to be exact in all dimensions [l7[ 
and is thus 'superuniversal'. For p < Pc{w), the mass is homogeneously distributed in the system with a mass 
distribution that has an exponential tail for large mass. At p = Pc{w), the mass distribution decays as a power law 
and for p > Pc, sl condensate forms on a single site that carries the additional macroscopic mass (p— Pc)L and coexists 
with a critical background fluid [l^ 

p(m) ~ exp[— m/m*] for p < pdw) 

~ ^ for p = pc{w) 
nv 

^ — — [-"condensate" for p > Pc (6) 

where the exponent r — 5/2 within the mean field theory [l| and is conjectured to have the same mean field value 
even in one dimension [l^ . 

Note that unlike ZRP, the exact joint distribution of masses P{mi,m2, ■ ■ ■ , wl) in the steady state is not known 
for the symmetric chipping model. In fact, it is believed [iTj that P{mi, 7712, . . • , 777^) does not have a simple product 
measure (factorisable) form as in ZRP in Eq. Another important difference is that in ZRP, the condensation 

transition happens both for asymmetric as well as symmetric transfer of masses to the neighbours, as long as one 
chooses the rate U{m) appropriately. In contrast, for the chipping model, a true condensation transition happens 
in the thermodynamic limit only for the symmetric transfer of masses. For the asymmetric transfer of masses (say 
only to the right neighbour) , the condensed phase disapears in the thermodynamic limit even though for finite L one 
does see a vestige of condensation transition [18]. However, a generalization that includes both the chipping model 
and ZRP as special cases does appear to have a condensation transition even with asymmetric hopping 19]. Finally, 
when the diffusion rate depends on the mass of the departure site in the chipping model in a certain manner, the 
condensation transition disappears [l^, [2l| . 

This simple chipping model with aggregation and fragmentation rules have been useful in various experimental 
contexts such as in the growth of palladium nanoparticles [l^. Besides, the possibility of such a condensation 
phase transition driven by the aggregation mechanism has been discussed in a system of Au sputtered by swift h eavy 
ions [23j . Finally, the chip ping model and its various generalizations have also been studied in the context of traffic [2J] , 
finance [23 and networks 12611. 



C. Asymmetric Random Average Process 

Another simple mass transport model that has been studied extensively [13, [H, [2^, [s^] is the asymmetric random 
average process (ARAP). As in the previous two models, ARAP is defined on a one dimensional lattice with periodic 
boundary conditions. However, in contrast to ZRP and the chipping model, here the mass 771^ at each site i is assumed 
to be a continuous positive variable. The model has been studied for both continuous-time as well as discrete-time 
dynamics [l^, [13] . In the continuous-time version, the microscopic evolution rules are [13, [11] : 

• In a small time interval dt, a random fraction 777,^ of the mass 771^ at site i is transported to the right neighbour 
{i + 1) with probability dt, where G [0, 1] is a random number chosen, independently for each site i, from a uniform 
distribution over [0, 1]. 

• With probability 1 — dt, nothing happens. 

The dynamics evidently conserves the total mass M — pL. At long times, the system reaches a steady state. Once 
again, the joint distribution of masses P(77ii, 7772, . . . , 777^) in the steady state does not have a factorised product 
measure form as in ZRP. 

How does the single site mass distribution p{m) look like in the large L limit? First important point one notices 
here is that the density p obviously sets the overall mass scale in this model. In other words, the mass distribution 
p{m, p) for any given p must have an exact scaling form 

1 / 777 \ 

p(777,p) = -Ff-j (7) 
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where the scaling function F(x) must satisfy conditions 

F{x)dx = l; I xF{x)dx = l. (8) 
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The first condition follows from normalization, J p{m, p)dm = 1 and the second from the mass conservation, 
J mp(ra, p) dm = p. Since the dynamics involves transfering a uniform fraction of mass from one site to its neighbour, 
the scaling in Eq. ([7]) is preserved by the dynamics. This is in contrast to ZRP or the chipping model, where one 
chips of a single unit of mass from one site to its neighbour and thereby the dynamics introduces a separate mass 
scale (unit mass), in addition to the overall density p. 

The scaling function F{x) for ARAP has been computed within the mean field theory [27l. [28l| 

F{x) = ^e-^/^ (9) 

and this mean field result is remarkably close to the numerical results in one dimension, even though one can prove 
rigorously 28] that the joint distribution of masses do not factorise. In contrast, for ARAP defined with a parallel 
discrete-time dynamics (where all sites are updated simultaneously), it has been proved [27l. [28| that the joint dis- 
tribution of masses factorise as in ZRP and the scaling function F{x) for the single site mass distribution can be 
computed exactly, F{x) = Axe^^^. The steady state of the discrete-time ARAP is also related to the steady state of 
the so called g-model of force fluctuations in granular materials [3l[ 

What about condensation? In ARAP, one does not find a condensation transition. This is of course expected since 
the density p just sets the mass scale and one does not expect to see a change of behavior in the mass distribution upon 
increasing p, apart from a trivial rescaling of the mass at all sites by a constant factor p. However, one can induce a 
condensation transition in ARAP by inducing an additional mass scale, e.g. , by imposing a maximum threshold on 
the amount of mass that may be transferred from a site to its neighbour [30j . 



III. A GENERALIZED MASS TRANSPORT MODEL 



Let us reflect for a moment what we have learnt so far from the three models discussed above. It is clear that the 
dynamics of mass transport often, though not always, may lead to a steady state that exhibits real-space condensation. 
For example, the ZRP and the symmetric chipping model exhibit real-space condensation, but not the ARAP. Also, 
we note that some of these models such as ZRP have a simple factorisable steady state. But factorisability of the 
steady state is clearly not a necessary condition for a system to exhibit real-space condensation, as we have learnt, 
e.g., from the study of the symmetric chipping model where the steady state is not factorisable. The factorisability, 
if present, of course helps the mathematical analysis. 

So, a natural question is: given a set of microscopic mass transport rules, what are the necessary and sufficient 
conditions that they may lead to a steady state that exhibits real-space condensation? For example, from the study 
of the three models above it seems that in order to have a condensation one needs to introduce via the dynamics a 
different mass scale, in addition to the density, such as the chipping of a single unit of mass in ZRP and the symmetric 
chipping model, or via introducing a maximum cap on the mass to be transferred in ARAP. If there is only one overall 
mass scale (density) that is preserved by the dynamics as in the usual ARAP, one does not expect to see a phase 
transition in the mass distribution as one changes the density. 

This question about finding the conditions for real-space condensation in a generic mass transport model seems too 
general and is perhaps difhcult to answer. Instead, one useful strategy is to restrict ourselves to a special class of mass 
transport models that have a factorisable steady state and then ask for the criterion, mechanism and nature of the real- 
space condensation phenomenon within this restricted class of mass transport models, which includes ZRP as a special 
class. This strategy has been demonstrated to work rather successfully in a recent series of papers [H, HI, [sj, [H, HI] 
and a fairly good understanding of the real-space condensation phcnonemon has been developed within this restricted 
class of mass transport models. This is what we will briefly discuss in this section. 

A Generalized Mass Transport Model: One can include all the three models discussed in Section II in a 
more generalized mass transport model [32] . For simplicity we define the model here on a one dimensional ring of L 
sites with asymmetric transfer rules, but it can be generalized in a straightforward manner on arbitrary graphs and 
arbitrary transfer rules. Similar mass transport models with open boundaries as well as dissipation at each site have 
also been studied [13] , though here we restrict ourselves to periodic boundary conditions and non-dissipative dynamics 
that preserve the total mass. On each site of the ring there is a scalar continuous mass mi. At any given time t one 
chooses a mass < Hi < mi independently at each site from a probability distribution (l>{pi\mi), normalized such that 
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JJ" (/)(/x|m) dfi = 1. In the time interval [t,t + dt] , the mass fXi is transferee! from site i to site i + simultaneously for 
all sites i (see Fig. ([1])). In a ring geometry, one identifies the site (L + 1) with the site 1. Thus, after this transfer, 
the new masses at times t + dt are given by [s^] 

m,(t + dt) = m,{t) - n,{t) + ^l^-l{t) (10) 

where the second term on the right hand side denotes the mass that has left site i and the third term denotes the mass 
that came to site i from site (i — 1). The function (j}{^\m) that specifies the distribution of the stochastic mass to be 




FIG. 1: Generalized mass transport model where a random mass /ii is transfered from site i with mass rrii to site i + 

transfered from any given site will be called the 'chipping kernel'. Here we take a homogeneous chipping kernel <j)[^\m) 
which does not depend on the site index. Note that the model above has been defined with parallel dynamics where 
all sites are updated simultaneously. Of course, by chossing dt — > 0, one can recover the continuous-time random 
sequential dynamics where the probability that two sites will be updated simultaneously is very small ^ 0{{dt)^). 
Thus the parallel dynamics includes the continuous-time (random sequential) dynamics as a special case. Note that 
for random sequential dynamics 0(/i|m) must generically be of the form 



4>{p\m) = a{^j,\'m)dt + 



1 — dt a{ii\ra)d^' 







5{p) (11) 



where a(/i|m) denotes the rate at which a mass /i leaves a site with mass m and the second term denotes the probability 
that no mass leaves the site. The form in Eq. pT|) is designed so that it automatically satisfies the normalization 
condition: (f){p\m) d/i = 1. 

This model with a general chipping kernel 0(/i|m) includes the previously discussed ZRP, chipping model and ARAP 
as special cases [32,] . Since we introduced these models in the previous section in continuous time, we will consider 
here the generalized model with chipping kernel of the form in Eq. (|lip with a general chipping rate a(/i|m). But of 
course one can consider a more general 4){fj,\m) with parallel dynamics that includes the continuous-time dynamics as 
a special case. Let us consider the three examples: 

(i) As a first example, we see that ZRP is recovered if in Eq. (fTT|) we choose for < ^ < m 

ain\m) = U{m)S{n-l) (12) 



Note that U{m) is zero if m < 1. 
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(ii) Similarly, the asymmetric chipping model is recovered is we choose 

a{fi\m) — wS{^ ^ 1) + S{fJ- ^ (13) 

where the first term refers to the event of a transfer of single unit of mass with rate w, the second term refers to the 
transfer of the full mass m with rate 1. 



(iii) Finally, one recovers ARAP by choosing 



a(/i|m) = — (14) 
m 



for all < /i- < m corresponding to the transfer of a fraction of mass that is chosen uniformly in [0, 1] leading to a 
uniform rate a(/x|m) independent of ^. 

By appropriately choosing the chipping kernel 0(/z|m), or equivalently the rate a{ii\m) for continuous-time dynamics, 
one can construct a whole class of mass transport models thus justifying the name 'generalized mass transport model'. 

Given a general chipping kernel 0(/i|m), or equivalently the chipping rate a{p,\m) for continuous-time dynamics, one 
can ask two important questions: (i) What is the steady state joint mass distribution P{mi, m2, . . . , tol)? (ii) Which 
types of (/)(/x|m), or equivalently a{p\m) for continuous-time dynamics, lead to a real-space condensation transition in 
the steady state? As discussed earlier, answers to either of these questions are hard to provide for a general chipping 
kernel 0(^|m) (or chipping rate a{^\m)). However, let us now restrict ourselves only to those chipping kernels 4>{^\m) 
(or a{^\m)) that lead to a factorised steady state distribution of the form 




P(mi,m2, . . . ,mi,) = ^^-^^/(mi) /(ma) . . . /(tol) (5 | } ^rrn-M] (15) 

where the partition function Zl{M) is just a normalization constant. This then leads us to the restricted mass 
transport model with factorisable steady state and we will address the questions regarding real-space condensation 
within this restricted class. The answers to these questions turn out to be easier for this restricted class since one can 
make use of the exact form of the steady state joint mass distribution (|15p . 



A Restricted Mass Transport Model: Here we restrict ourselves only to those chipping kernels that lead 
to a factorisable steady state (jl5p . Let us first investigate the question: given 0(/x|m) (or equivalently a{^\m) for 
continuous- time dynamics), what is the necessary and sufficient condition on (/)(/i|m) that leads to a factorisable 
steady state as in Eq. (jl5[) and if it happens, what is the exact form of the weight function /(m) in terms of (/)(/i|m)? 
Fortunately, answers to both questions can be obtained exactly which we state below without giving the details of 
the proof (see Ref. ([12]) for details). The necessary and sufficient condition for factorisability is that 0(/i|r7i) must 
be of the form [s^l 

nf^\m) = rm , . , --7- (16) 

Jo v[y)w[m-y)dy 

where v{x) and 'w{x) are arbitrary positive functions and the denominator is chosen to ensure normalizability, 
Jp (j){fi\'m)dfi = 1. In other words, if the chipping kernel 0(/i|m) factorises into a function of the mass that leaves the 
site and a function of the mass that stays on the site, then the steady state is guaranteed to be factorisable as in Eq. 
(lisp with a weight function whose exact form is given by the denominator in Eq. (|16p 



/(m) = / v{y)w{m-y)dy. (17) 
Jo 

This is the sufficiency condition. On the other hand, this condition can also be proved to be necessary, i.e., given that 
the steady state is of the form (fT5)l with some f{m), the chipping kernel has to be of the type (|16p and /(m) then 
must have the form ([T7]). 

Analogous condition can be found for continuous-time dynamics where (j){fj,\m) has the form pip with a chipping 
rate a(/x|m). The necessary and sufficient condition stated above for (/)(/z|m), translates into the following condition 
on a(/i|m) [3^ . 

a(M|m) = yi^i) ^ (18) 
z(m) 
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where y{x) and z(x) are two arbitrary positive functions. Note that a(fi\m) is a rate (and not a probabihty) and hence 
there is no normHzation condition here to be satisfied. If the rate a(/i|m) has the form (jlSp then we are gauranteed 
to reach a factorisable steady state (fTSl) with a simple weight function [34I 

/(m) = z{m). (19) 



As an example, it is easy to verify that the chipping rate in ZRP (|12p can indeed be written in the form in Eq. ([T| 
by choosing ?;(/i) = 5{ji— 1) and z[m) = HfcLi TjJ£) ^'-'^ — '^('^) ~ -'^ ^("^ < 0) = 0. Thus, ZRP with sequential 
dynamics is gauranteed to have a factorisable steady state (|15p with the weight function, /(m) = z{m) = Y[T=i TjIt) 
for m > 1 and /(O) = 2;(0) = 1. In contrast, both for asymmetric chipping model and the ARAP with chipping rates 
given respectively in (|13p and can not be written in the form (fTS]) with some choice of nonnegative functions y{x) 
and z{x), proving that neither of these two models has a factorisable steady state. 

One can ask several other interesting related questions. For example, suppose we are given a chipping kernel 
(/)(/x|to) and we want to know if it has a factorisable steady state or not. This amounts to an explicit search for 
suitable nonnegative functions v{x) and wlx) such that the kernel 0(/i|m) can be written as in Eq. (jl6p . This is often 
laborious. Can one devise a simple test which will allow us to do this search quickly just by looking at the functional 
form of (j){n\m)7 It turns out that indeed one can devise such a simple test which is stated as follows js^]: Given 
(/)(/x|m), first set m = fj, + a and compute the following two derivatives 

q{ti,a)^d,,d^log[cj){p\ti + a)] (20) 

which, in general, is a function of both variables fi and a. The test devised in Ref. states that a given 0(/^|m) will 
lead to a factosibale steady state iff this function q{fi, a) in Eq. ([^ is only a function of the single variable /U + cr, 
i.e., 

qi^i,a) =h{^i + a) (21) 
and in that case the weight function /(m) in the factorisable steady state (fT5|) is given explicitly by 



/(to) — exp 



dx / dyh{y) 



(22) 



In the discussion above, we have focused only on the ring geometry with asymmetric transfer of mass. Some of these 
results can be partially generalized to higher dimensions and even to arbitrary graphs [H, [11] and also to transport 
models with more than one species [39| of scalar variables, such as mass and energy for instance. 



IV. CONDENSATION IN MASS TRANSPORT MODELS WITH FACTORISABLE STEADY STATE 

In this section we discuss issues related to condensation within the restricted class of mass transport models that 
have a special steady state-namely a factorisable joint distribution (fTSjl with a suitable weight function f{m). We 
note that some aspects of the condensation transition in such a factorisable steady state, notably properties in the 
fluid state, was first studied in the context of the Backgammon model ^40i] without recourse to the dynamics that 
gives rise to such a factorisable steady state. A more complete analysis including the study of the properties of the 
condensed phase was undertaken in Refs. (3^ . [s^ which will be summarized in this section. 

There are three main issues here: (i) criterion: what kind of weight functions f{m) lead to a condensation transition 
(ii) mechanism: what is the mechanism of the condensation transition when there is one and (iii) nature: what is 
the nature of the condensate, e.g., what is the statistics of the mass in the condensate in the condensed phase? All 
these questions can be answered in detail for factorisable steady states. We briefly mention the main results here, the 
details can be found in Rcfs. [s^. [s^. 

Crietrion: The factorisation property allows one to find the criterion for a condensation transition rather easily 
by working in the grand canonical ensemble (GCE). Within GCE framework, one introduces a fugacity exp[— sm] 
where s is the negative of the chemical potential associated with each site. This is just equivalent to taking the Laplace 
transform of Eq. ((TSl) with respect to the total mass M (with s being the Laplace variable), which replaces the delta 
function by exp[— s(toi + TO2 + . . .to^)]. Then s is chosen such that the constraint M = ^ is satisfied on an 
average. Given that each site now has a mass distribution p(to) — /(TO)exp[— sm] (upto a normalization constant), 
the equation that fixes the value of s for a given M — pL is simply 
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The criterion for condensation can be derived easily by analysing the function p(s) defined in Eq. (|23|) . If for a 
given p, one finds a solution to this equation s = s* such that the single site mass distribution is normalizable, i.e., 
/ p{ni) dm = J /(m) exp[— s* m] dm is finite, then there is no condensation in the sense that for all values of p, the 
single site mass distribution has an exponential tail and there is not one special site that needs to accomodate extra 
mass. On the other hand, it may be that for certain /(r7i)'s, as one increases p, there may be a critical value pc below 
which one finds a good solution s to Eq. (|23|) . but such a solution ceases to exist for p > pc- This will then signal the 
onset of a condensation because for p > pc, the system needs to break up into two parts: (a) a critical background 
fluid part consisting of (L — 1) sites at each of which the average density is critical pc and (b) a single condensate site 
which accomodates the additional mass {p — Pc)L. 

As an example, let us consider f{m) that decays slower than an exponential, but faster than l/m^ for large m. 
Since f{m) decays slower than an exponential, in order that the single site mass distribution p{m) — f{m) ™ is 
normalizable (i.e., / p{m) dm = 1), the possible solution s* of Eq. (|23p can not be negative. Thus the lowest possible 
solution is s* — 0. Now as s ^ 0, the function p{s) in (|23p approaches a critical value. 



Pc = p{s 0) 



/o°° ™ /(m) dm 



(24) 



which is finite since /(m) decays faster than for large m. Thus as long as p < pc, by solving ((23)) one will get a 

positive solution s* and hence no condensation. As p — > pc from below, s* — > from above. But for p > pc, there is 
no positive solution s* to (j23p . which signals the onset of a condensation transition. 

A detailed analysis of Eq. shows [35] that in order to have condensation, the weight function /(m) must have 
a large m tail that lies above an exponential but below i.e., exp[— cm] < /(m) < for large m with some 

positive constant c > 0. A natural candidate satisfying this criterion is 



/(m) ~ A', 



with 7 > 2 



for large m. Indeed, the ZRP discussed in the previous sections with the choice U{m) ^ (1 
to a weight function f{m) \2b\ and then condensation happens only for 7 > 2. 



(25) 

7/to) for large m leads 



Mechanism and Nature: Given an appropriate weight function /(m) such as in ([^5]) that leads to condensation, 
one can then ask about the mathematical mechanism that drives the condensation. Actually, there is a very simple 
way to understand this mechanism in terms of sums of random variables [3^] which we will discuss in the next section. 
For now, we notice that in an infinite system, where GCE is appropriate, the single site mass distribution p{m) has 
the form p{m) = f{m) exp[— sm] with an appropriate s which is the solution of Eq. (|23p as long as p < pc. For p > pc, 
there is no solution to Eq. ([23l) . In fact, for p > pc, the value of s sticks to its critical value Sc and the GCE framework 
is no longer valid. To understand how the condensation manifests itself in the single site mass distribution, one has 
to study a system with a finite size L and work in the canonical ensemble with a strict delta function constraint as in 
Eq. (USD. 

In a finite system of size L, the single site mass distribution p{m) can be obtained by integrating the joint distribution 
(fT^ over the masses at all sites except one where the mass is fixed at m. It is easy to see from Eq. p5|) that 



pirn) 



P(m, mi, m2, . . . mi)dm2 dm^ . . . dm^ = f{m 



ZL-i{M-m) 
Zl (m) 



where the partition function Zl{M) 



Zl{M) ^ J f{mi)f{m2).f{mL)5 {^Yl - c'™! dm2 ■ ■ ■ dmL- 



Taking Laplace transform with respect to M, one gets 

ZL{M)e 



dM = 



/(m)e^"™ dm 



(26) 



(27) 



(28) 



which can be formally inverted using the Bromwich formula 



Zl{M) 



27r7 



exp [L {lng{s) + ps)] 



(29) 
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where we have used M = pL, the integral runs along the imaginary axis (sq + * y) in the complex s plane to the right 
of all singularities of the integrand and 

poo 

g{s)= / /(m)e-^"dm. (30) 
Jo 

One can write a similar integral representation of the numerator Zl(M — m) in Eq. (|26p. Next one analyses Zl{M) 
and Zl{M — m) using the method of steepest decent in the large L limit. 

As long as there is a saddle point solution to Eq. ((29|l . say at s = s*, we see immediately from Eq. ([26l) that the 
single site mass distribution for large L has the form, p{m) ~ /(m) exp[— s*m], i.e., one recovers the GCE result. 
Thus the GCE aprroach is valid as long as there is a saddle point solution s* . As one increases p, the saddle s* starts 
moving towards in the complex s plane and when p hits p^-, s* 0. For p > p^, there is no saddle point and one has 
to analyse the Bromwich integrals by correctly choosing the contour to evaluate p(m) for p > pc- This was done in 
details in Ref . [3^] . We omit the details here and mention the main results for Zl {M) in Eq. (P^j) and subsequently 
for p{m) in Eq. ((26|) . when the weight function /(m) is chosen to be of the power law form ((25)l . 

If we normalize f{m) such that f{m)dm = 1, the partition function in Eq. (|29p can be interpreted as the 
probability that the sum of L i.i.d variables, each drawn from a distribution f{m), is M (see Section V also). 
Analysing the Bromwich integral (|29p for large L, using the small s behavior of g{s) in Eq. ((301) . on<5 finds that 
the asymptotic behavior of the distribution Zl{M) is different for 2 < 7 < 3 and for 7 > 3. 

2 < 7 < 3: In this regime, one finds the following scaling behavior of Zl{M) 

'p^L~M' 



where Pc — fJ-i — /q°° m f{m) dm is the first moment and the function V^{z) is given explicitly by [3£ 



(31) 



V-f{z) = - / dy e-''^y''~\oa [bcos{T:-f /2)y'<~'^ + yz] . (32) 
Jo 

Here C3 — — 6sin(7r7/2) > and b — ^r(l — 7) for 2 < 7 < 3 with A being the amplitude in ((25)) . The precise 
asymptotic tails of this scaling function can be computed 35] 

V^{z) ~ A|zr^ asz-^-00 (33) 
= Co at z = (34) 

^ Ciz(3-7)/2(7-2)g-c..<-^)/<-^) aSZ^OO (35) 

where cq, ci and C2 are known constants 35]. Thus the function is manifestly non-gaussian. 
7 > 3: In this regime, the partition function Zl{M) has a gaussian peak 

Zl{M) ~ ^ g-(M-p.L)V2A^L _ ^^^1 ^ o{L^'^) (36) 

V 27rA2 L 

where ~ p2 ~ /i? with fik = /q°° "m^ /(™) dm being the /c-th moment. But far to the left of the peak, Zl{M) has 
a power law decay [35| . 

So, how does the single site mass distribution p{m) in Eq. ((26]) look like? We have to use the result for the partition 
function derived above in Eq. ((26|) . We find different behaviors of p{m) in different regions of the (p — 7) plane. For 
7 > 2, there is a critical curve pdj) in the (p — 7) plane that separates a fluid phase (for p < pdj)) from a condensed 
phase (for p > Pc{l))- In the fluid phase the mass distribution decays exponentially for large m, p{m) ~ exp[— m/m*] 
where the characteristic mass m* increases with increasing density and diverges as the density approaches its critical 
value Pc from below. At p = pc the distribution decays as a power law, p{m) ~ m~^ for large m. For p > pc, the 
distribution, in addition to the power law decaying part, develops an additional bump, representing the condensate, 
centred around the "excess" mass: 

Me, = M-pcL. (37) 



Furthermore, by our analysis within the canonical ensemble, we show that even inside the condensed phase (p > Pc(7)), 
there are two types of behaviors of the condensate depending on the value of 7. For 2 < 7 < 3, the condensate is 




FIG. 2: Schematic phase diagram in the p-7 plane. 




FIG. 3: The exact single-site mass distribution p{m) plotted for a particular choice of /(m) with 7 = 5/2 and pc = 2, and 
system size L = 100: full line p = 1 (subcritical: fluid phase); dotted line p — 2 (critical); sadhes line p = 6 (supercritical: 
condensed phase), the condensate bump, Pcond{ni), is evident in the supercritical phase. 
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characterized by anomalous non-gaussian fluctuations whereas for 7 > 3, the condensate has gaussian fluctuations. 
This leads to a rich phase diagram in the (p — 7) plane, a schematic picture of which is presented in Fig. 

Detailed form of p{m) for p < Pc (fluid phase) , p > pc (condensed phase) and p = Pc (critical point) are summarized 
below (see Fig. ^ also). 

Fluid phase p < p^- In this case one finds 

p(m) - /(m)e-™/'"' for l<m<M (38) 

where the characteristic mass m* diverges p approaches pc from below as (p — Pc)~^ foi' 7 > 3 and {p ~ Pc)"^^^'^^'^^ 
for 2 < 7 < 3. 

Condensed Phase p > p^. In this case one finds 

p{m) ~ /(to) for K TO < 0{L) (39) 

p{m) ~ f{rn)- — for to = xAI^x where < x < 1 (40) 

(1 — xy 

p{m) ^ Pcond(™) for TO (41) 

Here Pcond is the piece oi p{m) which describes the condensate bump (see Fig. ([3]): Centred on the excess Mex and 
with integral being equal to l/L, it takes on two distinct forms according to whether 7 < 3 or 7 > 3. 
For 2 < 7 < 3 



Pcond 



(to) ~ 



Li/(7-i) 



(42) 



where the function V-y{z) is given explicitly in (|32p . Thus clearly the condensate bump has a non-gaussian shape for 
2 < 7 < 3 and we refer to this as an 'anomalous' condensate. 
On the other hand, for 7 > 3 

V2ttA^L^ 

i-e. Pcond(»7i) is gaussian on the scale \m—Mex\ ^ 0{L^^^), but, far to the left of the peak, p{m) decays as a power law. 
Critical density p ^ pc- In this case one finds that 

p{m) cx fim)Vy (m/L^/^''-'^^^ , for 2 < 7 < 3 (44) 

p{m) cx /(to) e-'"'/^^'^ 7>3. (45) 

where the scaling function V-y{z) is as before. Thus at criticality p(rn) decays as a power law m~'^ for large to which 
is cut-off by a finite size scaling function and the cut-off mass scales as 

mcut-off' - for 2 < 7 < 3 (46) 

~ foj. ^ > 3^ (47) 

Physical picture: It is useful to summarize the main physical picture that emerges out of this mathematical analysis. 
We notice from the joint distribution (jlSp that the masses at each site are 'almost' independent random variables 
each with a power law distribution /(to), except for the global constraint of mass conservation imposed by the delta 
function which actually makes them 'correlated'. The system feels this correlation for p < pc and p > pc m different 
ways and exactly at the critical point p — pc the effect of the constraint is actually the least. For p < pc^ the effective 
mass distribution at each site acquires an exponential tail, p{m) ^ f{m) exp[— s*m] which is induced by the constraint. 
For p = Pc s* — > and p{m) ~ f{iTi), thus the system does not feel the constraint at all and the masses behave as 
completely independent random variables each distributed via /(m). But for p > pc, while (L — 1) sites behave as the 
critical fluid, i.e., as if the mass at each of these (L — 1) sites is distributed via /(to), there is one single condensate 
site which acquires the additional mass (p — Pc)L. 

For p > Pc, the resulting non-monotonous shape of the single site mass distribution (with an additional bump) in 
Fig. ^ can then be understood very easily from this physical picture. Basically, for p > pc, the total mass M of the 
system splits into the critical fluid and the condensate part: 

A^fiuid (48) 
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where mcond denotes the mass at the condensate and the critical background fluid mass 

L-l 
i=l 

is a sum of [L — l) independent random variables (masses) each distributed via f{m). Thus the probability distribution 
-^^fluid is given precisely by the partition function Z^^i (Mfluid) in Eq. (j27p . This partition function can be computed 
explicitly and the results are given in Eqs. (I31|) and ()36l) . Knowing this partition function, the distribution of the 
condensate mass mcond can be obtained using Eq. (|48p giving 

Prob(mcond = y) = ^L-i(Af-y). (50) 

The overall single site mass distribution for p > then can be computed as follows: if we choose a site at random, 
with probability {L — i)/L it belongs to the background fluid and hence its mass distribution is f{m) whereas with 
probability l/L it will be the condensate site with mass distribution given in Eq. (|50p . Thus for p > pc 

p(jn)^^^^j^f{m) + jZL-i{M-m) (51) 

The second term is what we referred to before as 

Pcond(m) = ZL-iiM - to) (52) 

and it is this piece that describes the bump in Fig. Q for p > pc- Its precise asymptotic behavior is detailed in Eqs. 
dUl) and ^ respectively for 2 < 7 < 3 and 7 > 3. 



V. INTEPRETATION AS SUMS AND EXTREMES OF RANDOM VARIABLES 

There is a very nice and simple way [lH, using sums of random variables, to understand the mechanism of the 
condensation transition for factorsable steady states (fTS)) with a given weight function /(to) say of the form (|25p . 
Let us consider a set of L positive i.i.d random variables {toi, TO2, . . . , tol} each drawn from a distribution /(to) (we 
choose /(to) such that it is normalized to unity). Let M = X]i=i "^i be the sum. For instance, M can be interpreted 
as the position of a random walker after L independent steps of lengths mi, m2, . . ., m^. Then we notice that the 
partition function Zl{M) in Eq. (^7)) can be interpreted as the probability that the walker reaches M in L steps 
starting from the origin. 

How does one interpret condensation in this random walk language? Note from Eq. (j24p that for normlaized (to 
unity) /(m), the critical density 

/•oo 

Pc = Ml = / mf{'m)dm (53) 
Jq 

is just the mean step length of the random walker's steps. Thus if the final position M < pcL = piL (i.e., p — M/L < 
Pc), we would expect that the typical configuration of the random walker's path would consist of steps each of which 
is of ^ 0(1)- But, for AI > piL (i.e., p > pc), the ensemble will be dominated by configurations where {L ~ 1) steps 
are of ~ 0(1), but one single big step of ~ (M — pi{L — 1)) ^ {p — Pc)L to compensate for the deficit distance. 
This single big step is precisely the condensate. Within this interpretation, it also becomes clear that for /(m) of 
the form (|25p . there are two possibilities depending on whether the second moment of the step length distribution 
P2 = /q°° fiiTi) dm is divergent (2 < 7 < 3) or finite (7 > 3). In the former case, the corresponding random walker 
is a Levy flight with analomalosly large fluctuations that leads to an anomalously large fluctuation in the condensate 
mass. In the latter case, by virtue of the central limit theorem, one recovers a gaussian fluctuation leading to a 
gaussian distribution of the condensate mass. This explains the two types of condensate phases in the phase diagram 
in Fig. (0). 

Condensation and Extreme Statistics: Another interesting issue intimately related to the condensation is the 
associated extreme value statistics-e.g, what is the distribution of the largest mass TOmax in the system, where 

mmax = max(mi,m2, . . . ,tol). (54) 

This is particularly important in the condensed phase where the largest mass is carried by the condensate, at least in 
models with simple factorised steady states where there is a single condensate. The theory of extreme value statistics 
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is well developed in cases where one studies the extreme (e.g., the maximum) of a set of i.i.d random variables [4l|. 
In our case, the factorised steady state in Eq. (|15p shows that the masses are not completely independent, but are 
correlated via the global mass conservation constraint explicitly manifest in the delta function in Eq. (fT5| . Without 
this delta function constraint, with /(m) m~'^ and 7 > 2, the scaled distribution of the maximal mass would have 
been Frechet distribution [4l|. However the presence of the constraint induces important correlations that changes 
the nature of the distribution of the maximal mass. 

In Ref. [4^, the authors studied rigorously how the typical value of the extremal mass scales with system size 
in ZRP with 7 > 3. For more general mass transport models with factorised steady states as in Eq. (jlSp. the full 
distribution of extremal mass was studied recently in Ref. j43|. In was found that in the fluid phase (p < pc), 
where the single site masses effectively become uncorrelated but with an additional exp[— s*m] factor that comes from 
the conservation constraint: p{m) ~ /(m) exp[— s*to] and hence the maximal mass, in the scaling limit, becomes 
Gumbel [i^. At the critical point where s* = (i.e., where the constraint is least effective), one recovers the Frechet 
distribution. But for p > pc the maximal mass distribution is the same as that of the condensate mass, i.e.. 

However, as mentioned in the previous section, the distribution of mcond = M — Mfiuid can be computed via computing 
the distribution of Mfluid as a sum of independent random variables each distributed via f{m) and is given by Eq. 
([50|. Thus, 

Prob (m,„ax = y) ~ Zl-i{M - y) = Lpcondiv) (56) 

where Pcond 

(m) is given in Eqs. (02]) and (|43]) respectively for 2 < 7 < 3 and 7 > 3. Thus, in the condensed 
phase, one has a completely new type of extreme value distribution of correlated random variables which is exactly 
computable 'is']. Moreover, it is interesting to note that for p > pc, the computation of the distribution of the extreme 
of correlated variables reduces to the calculation of the distribution of sum of independent random variables. 



VI. CONCLUSION 



In this brief review I discussed recent developments in understanding the physics of the real-space condensation in 
a class of mass transport models. Here, real-space condensation refers to the phenomenon when, upon increasing the 
density beyond a critical value, a macroscopically large mass settles onto a single site in real space in the steady state. 
The system discussed is homogeneous in the sense that the transport rules do not depend on the sites or particles. 
Thus, in the limit of an infinite system size, the formation of the condensate at one single site actually breaks the 
translational symmetry in the system spontaneously. The criterion and mechanism of the transition as well as the 
detailed finite size dependence of the distribution of the condensate mass in the condensed phase was discussed within 
a restricted class of one dimensional mass transport models that have a factorisable steady state. 

There are several directions in which the questions addressed here can be extended, some of which are briefly 
mentioned below. 

Pair factorised steady states: Here we discussed only mass transport models that have a factorised steady state 
(lisp . A natural question is whether real-space condensation can happen in other types of steady states that are not 
simply factorisable as in (jl5p and if so, (i) are there natural local transport rules that lead to such steady states and 
(ii) does the nature of the condensate change fundamentally from the one with factorisable steady states? Recently, 
a generalization of Eq. (fTS]) . called the pair factorised steady states (PFSS) was introduced in Ref. [3] 

1 ^ / L \ 

F(mi,m2,...,TOL)] = -^--^^^5(m^,mi+i)(5 ( ^mi-Mj . (57) 

Thus there is one factor g{mi,mi^i) for each pair of neighbouring sites on a rin g of L sites. The transport rules, 
involving three neighbouring sites, that lead to such steady states were also found [4J]. Interestingly, the condensate 
in this PFSS, for a class of weight functions g{m, n) that are short ranged, was found to spread over a relatively large 
number of sites ~ 0(L^/^) \4^, in contrast to the condensate that forms over a single site in the usual factorised 
steady state P^ . The average shape of this sub-extensive condensate for a class of weight functions g{m,n) as well 
as the precise form of the 'condensate bump' in the single site mass distribution in the condensed phase were recently 
computed in a very nice paper [4^ . In addition, the transport rules that lead to PFSS on an arbitrry graph were also 
found recently [46], thus generalizing the one dimensional models with PFSS. 
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Dynamics: Here we have discussed only static properties associated with the condensation transition. Another 
interesting issue is the nature of the dynamics in the steady state as well as in the approach to the steady state, 
in particular in the condensed phase [23, H, [4^, [sO, [sij. In a finite system, the condensate forms at a site, then 
survives there over a long time T^, then desolves and then forms at another site. For ZRP with f{m) ~ with 
7 > 2 and in the stationary state, it was found that while the condensate life-time ~ {p — Pc)^"*"^ L'' for large 
L [4§|, there is another shorter time scale associated with the relocation of the condensate, Tr ^ (p — PcYL^ j50| . 
In addition, the current fluctuations in the steady state show a striking change of behavior [50] as one goes from the 
fluid phase {p < pc) to the condensed side (p > pc). On the fluid side, the current fluctuations shows an interesting 
oscillating behavior due to the presence of kinematic waves [53| , similar to the oscillatory behavior of the variance 
in the displacement of a tagged particle in 1-d asymmetric exclusion processs [s^ . Another interesting dynamical 
quantity is the power spectra associated with the time series depicting the evolution of the total number of particles 
over a fixed segment of a ring, studied recently in the context of ZRP [s^ . 

It would be interesting to study the dynamics for the mass transport models in higher dimensions or arbitrary 
graphs, and also for more generalized steady states such as the PFSS. 

Other interesting directions involve studying the condensation phenomena in multispecies models [isl . [39 L misan- 
thropic process fl3] , instability of the condensed phase due to nonconserving rates in the chipping model ll, IH^I as 
well as in ZRP[55i, ^5^,] and also in ZRP due to quench disordered particle transfer rate [s^l , ZRP leading to multiple 
condensates (5^, condensation in polydisperse hard spheres [11] etc. Most of these generalizations have been carried 
out so far in the context of ZRP, but it would be interesting to study the general mass transport model discussed 
here with these additional generalizations. It would also be interesting to compute the distribution of maximal mass 
in PFSS where the condensate is sub-extensive. In this case, the extremal mass is not the total mass carried by the 
full condensate, but rather the site inside the condensate that carries the largest mass. 
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